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Abstract 


The  diffraction  by  a  grating  is  examined  (for  spacing  large 
compared  to  wavelength  and  dimension  of  grating  elsment)  for  wavelengths 
in  the  neighborhood  of  the  "  Rayleigh"  wavelengths.  The  shape  of  the 
elements,  and  their  size  in  wavelengths  is  unrestricted.  The  results, 
including  the  effect  of  interaction,  are  expressed  in  terms  of  quanti- 
ties relating  to  single  scattering.  Some  properties  of  certain  deter- 
minants formed  from  single  scattered  amplitudes  are  derived.  The  re- 
sults are  compared  with  those  obtained  by  other  authors,  using  various 
restrictions  on  the  parameters. 


Table  of  Contents 

Pa^ 

1»  Introduction  1 

2,  Preliminary  results  3 

3«  The  case  of  resonance  6 

[(•  Remarks  on  energy  conservation  13 

5»  Comparison  with  other  theories  II4 

Appendices  22 

References  IjO 


-  1  - 


1.  Introduction 

The  purpose  of  the  present  work  is  to  present  the  results  of  an  in- 
vestigation of  the  behavior  of  diffraction  gratings  of  cylinders  for  wavelengths 
neighboring  certain  critical  wavelengths,  the  so-called  "  Rayleigh  wavelengths"* 
The  analysis  is  based  on  the  res-ults  of  a  prevloxis  investigation'--'**--',  in 
which  expressions  were  derived  for  the  Fourier  amplitudes  of  the  diffracted 
field  of  an  grating  of  arbitrary-  cylinders,  subject  to  the  restrictions  that  the 
spacing  is  large  as  compared  to  cylinder  dimension  and  wavelength,  and  that  the 
Rayleigh  wavelengths  are  excluded.  The  expressions  there  obtained  for  the 
Fourier  amplitudes  only  involved  quantitites  associated  with  a  single  scatterer, 
in  spite  of  the  fact  that  interactions  were  taken  into  account.  The  principal 
result  of  the  present  paper  is  the  extension  of  the  conclusions  of  [2]   to  the 
case  of  critical  wavelengths,  in  which  the  significance  of  interaction  is  maxi- 
mal. 

The  required  results  of  [l]  and  [2]   are  recapitulated  in  Section  2  of 
the  present  work.  They  form  a  special  case  of  a  general  theory^-^   for  the  cal- 
culation of  the  field  scattered  by  a  combination  of  obstacles  in  terms  of  the 
scattering  patterns  of  the  obstacles  in  isolation.  Experience  with  the  case  of 
diffraction  by  a  pair  of  half-planes  (i.e.,  a  slit)  '•  -^  indicates  that  the  ac- 
cxiracy  of  this  kind  of  calcvilation  is  good  down  to  spacing  of  one  wavelength. 

In  Section  3  of  the  present  paper  we  discuss  the  behavior  of  the  Fourier 
amplittide  of  the  grating  field  in  the  limit  in  which  the  critical  or  Rayleigh 
wavelengths  are  approached.  In  this  limit  the  "  single  scattering"  result  di- 
verges, while  the  result  given  in  Section  2  is  indeterminate.  However,  the  limit 
is  evaluated  explicitly  in  the  present  work.  We  also  include  in  Section  3  a  dis- 
cussion of  certain  determinants  involving  the  scattering  amplitudes  of  a  single 
cylinder;  the  behavior  of  these  determinants  is  crucial  for  the  limiting  process 
in  question. 


Section  I4  is  concerned  with  a  brief  discussion  of  energy  conservation. 
It  is  shovm  that  the  "  single  scattering"  result  conserves  energy  with  an  error 
which,  while  small  for  general  large  spacings,  contains  a  factor  which  diverges 
for  critical  spacings.  On  the  other  hand,  o\ir  results  are  free  of  this  defect. 

In  Section  5  we  compare  our  results  with  those  of  pre'/iotis  investigators. 
Most  previous  work  has  been  based  on  other  types  of  restrictions  than  qmt   own. 
However,  Twersty'--''   ,  who  treated  the  case  of  circular  cylinders  for  non- 
critical  spacings,  obtained  equivalent  results  for  the  contributions  of  the  in- 
dividual cylinders  in  that  case,  for  certain  angles  of  incidence;  and  these  re- 
sults were  extended  by  Wagnus  (in  Appendix  to  [5])  to  arbitrary  angles  of  in- 
cidence subsequent  to  the  developments  in  [l]. 

A  more  interesting  ccmparison  is  made  possible  by  considering  the  appli- 
cation of  our  results  to  the  reflection  grating,  which  is  the  seat  of  the  celebrated 
"Wood's  Anomalies"  discussed  by  Wood L' J,  RayleighLJ,  FanoL^-J  and  Strong'-  -'• 
A  definitive  theory  has  been  given  by  Artraann*-  -•.  Since  the  only  restriction  of 
Artmann's  analysis  is  that  the  corrugations  be  small  compared  to  wavelength,  a 
comparison  with  otcr  theory  is  possible  in  the  case  in  which  corrugations  are 
small  and  the  spacing  is  large.  We  specialize  our  results  to  small  corrugations 
in  Section  5,  and,  after  giving  explicit  approximate  formulas  for  the  scattering 
amplitudes  of  a  thin  strip  in  the  case  when  the  magnetic  vector  is  parallel  to 
the  generators  of  the  grating,  we  are  able  to  show  agreement  with  the  results  of 
Artmann. 

On  the  other  hand  it  is  pointed  out  that  both  for  the  resonance  and  non- 
resonance  cases  our  results  are  in  disagreement  with  results  given  in  a  new  re- 
port  by  Twersky^  -'.   The  nature  of  the  discrepancy  is  explained  and  possibili- 
ties for  its  removal  are  suggested. 

Certain  calculations,  omitted  in  the  main  body  of  the  text  in  order  to 
keep  the  argument  concise,  are  carried  out  in  detail  in  the  appendices. 
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2»  Preliininary  results 

It  is  well  known  that  in  the  far  field  every  outgoing  wave  function  takes 
the  form 
(1)         H^(kr)  f(e)  =  Y^   expCikr  -  i  J)  •  f(e)   , 

where  r  and  6  are  polar  coordinates.  Also,  in  a  problem  dealing  with  diffraction 
by  a  combination  of  obstacles.  Green's  theorem  shows  that  we  can  regard  the  scat- 
tered field  as  ccsnposed  of  a  sum  of  outgoing  wave  functions,  each  regtilar  in  the 
exterior  of  one  of  the  scattering  bodies.  Suppose  the  bodies  in  question  are 
sufficiently  far  apart  for  all  dimensions  of  each  body,  as  well  as  the  wavelength, 
to  be  small  compared  to  the  distance  to  any  other  body  in  the  array.  Then  it 
can  be  shown  from  (l)  that  in  the  neightborhood  of  any  body  in  the  configuration 
the  total  field  is  approximately  equal  to  the  incident  field  plxis  the  field  scat- 
tered by  the  body  in  question  plus  a  s\im  of  plane  waves  arriving  from  the  direc- 
tions of  the  other  bodies.   (The  complex  sunplitudes  of  the  latter  plane  waves 
are  at  first  unknown.]  Consequently  the  actual  field  contributed  by  each  body 
is  approximately  an  unknown  linear  combination  of  the  fields  which  that  body 
would  scatter  when  responding  to  the  plane  waves  ih  question  in  isolation.  The 
necessary  coefficients  can  then  be  calculated  by  a  consistency  argument.  For 
details  in  the  case  of  E||,  see  [2],  \}\   (the  argument  is  easily  seen  to  apply 
for  arbitrary  boundary  conditions,  however). 

The  theory  in  question  can  be  applied  to  an  infinite  grating  of  cylinders. 
Before  stating  the  results  we  recall  some  formulas  from  the  general  theory  of  in- 
finite gratings. 

We  consider  a  grating  of  congnient  cylinders  lying  in  the  plane  y  =  0» 
If  a  plane  wave  is  incident  at  the  angle  j2f,  so  that 

This  is  the  two-dimensional  resvilt.  Time  dependence  exp(-i(ot)  is  assumed. 


-  li  - 

(2)  u^  -  exp[ik(x  cos  ^  +  7  sin  ^], 

and  if  the  distance  between  corresponding  points  of  the  grating  is  d,  then  the 
fields  contributed  hj  the  various  cylinders  differ  among  each  other  only  by  a 
translation  and  by  a  phase  factor  expllknd  cos  ^j,   so  that  the  total  field  may 

be  written 

00 

(3)  u     -    u^^  +  21]    exp[ilmd  cos  0J  F(x-nd,y)   • 

-00 

Suppose  further  that  the  far  field  form  of  F(x,y)  is  given  by 

(U)        F(x,y)  -  /^  e:q?(ikr-i  J  )  f(©)  , 

where  x  «"  r  cos  ©,  y  ■  r  sin  Q»     Then  in  consequence  of  the  regularity  of  the 
field  in  the  exterior  of  a  horizontal  strip  containing  the  grating  elements, 
the  total  field  can  be  written  as  follows  (see  [2J): 

00 

-1 

■  00    -    "V 


(^)       ^  =  ^inc  -^  Z]   kdsine   •  ^^i  "v^   ^^[^^  ^°^^«  ^  ^v^l 

-OD  V  *-  -• 


except  for  values  of  y  in  horizontal  strips  containing  the  grating  elements.  The 
upper  sign  is  to  be  taken  when  y  is  aboTre  the  strip.  In  tliis  exact  formula  the 
spectral  angles  9    are  given  by  the  formula 

k  cos  6  »  k  cos  0  +   2nv/d  , 


k  sin  ©^  »  /k^-  (k  cos  0  +  2nv/d)' 


The  infinite  grating  problem  is  thus  eqxiiveilent  to  the  problem  of  de- 
tenaining  the  values  of  the  function  f(©)  for  ©  «  ©  .  A  finite  number  of  the 
quantities  ©  are  real;  these  lead  to  propagating  waves  in  (5)  (i.e.,  orders 
of  the  8pectnjm)o  The  function  f(©)  can  be  expressed  in  terms  of  the  boundary 
values  of  F(x,y)  by  quadratures  which  are  meaningful  also  for  ccanplex  ©.  We 
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can  now  state  o\ir  preliminary  formulas,  with  trtiich  we  operate  in  the  balance  of 
this  paper.  We  have  found,  using  the  type  of  analysis  described  in  the  foregoing 
paragraph,  that 
(7)  f(e)   -  fj(e)  +  Aj^f(0)fj^(9)  +  Aj^f(n)fj^(6)  , 

(7a)        r(x,y)  -  Fj(x,y)  +  A^fiO)Tj{x,y)   +  Ajjf(n)Fj^(x,y)   , 

where  the  fiinctions  f-j-(O),  ^t(®)  ^"d  foC©)  are  those  describing  the  far  field  of 
a  single  cylinder  in  the  grating,  when  excited  respectively  by  tiie  incident  field, 
a  plane  wave  coming  from  x  =  -oo,  and  a  plane  wave  coming  from  x  ■  +cd»  The 
fiinctions  F_(x,y),  etc.,  are  the  corresponding  scattered  fields  defined  for  all 
X  and  y»  The  constants  A^  and  A~   are  given  by  the  formulas 

00 


h     "      /S    ®^[-^f]    S    exp[inM  (1-cos  JJO]  / y^     , 
h    '      /S    ^["^f]   ^    exp[inkd  (1+cos  JJQ]  / y^    . 


(8) 

00 


The  quantities  f(0),  and  f(n)  (which  are  obtained  from  (7)  by  setting  6=0  and 
e  =  n  respectively,  and  sol'/ing  the  resulting  pair  of  simultaneous  equations), 
are  given  by  the  following  formulas 

1-1 


(9) 


f(0)  -  [f j(0)  -  A^a^]  |l  .  [A^fj^(O)  .  Aj^fj^(n)]  .  A^^A^^  aX 
f(n)  -   [f,(n)  -  A,cJ  |,  .  ^,^,^(,)  ,  ^,^(„)J  ,  ^,^  q^     ^ 


in  which  a. ,  Op,  a-  are  given  by  the  formulas 
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(10)  a^ 


fj(0)   fj^(O) 
fj(n)   fj^(n) 


fj^(C)   fj(C) 
fj^(n)   fj(n) 


f^iO)       fj^(n) 
fj^(O)   i^(n) 


This  result,  which  is  obtained  in  [2J,  holds  provided  that  k,  d,  and  ^  are  such 
that  there  is  no  positive  or  negative  integer  v  for  which  sin  ©  vanishes.  This 
implies  that  both  A,  and  A_  converge.  Then  A,  and  A^  are  small  qviantities  for 
sufficiently  large  kd,  so  that  the  leading  term  of  (7)  is  fj(©).  The  approxima- 
tion f(6)'^f-.(©)  thus  expresses  the  natural  assumption  that  each  grating  element 
responds  as  if  in  isolation,  when  the  spacing  is  very  large. 


3«  The  case  of  resonance 

Our  task  in  'che   present  work  is  to  discuss  the  behavior  of  the  solution 
in  the  neighborhood  of  those  frequencies  or  angles  of  incidence  that  mafce 
sin  6  =  0.  We  shall  refer  to  such  a  situation  as  a  "  resonance",  and  shall 
discuss  the  limiting  value  of  the  result  (7),  as  well  as  the  behavior  in  the 
neighborhood  of  the  limit,   ^e  shall  not  discuss  the  mathematical  question  of 
the  existence  and  properties  of  the  solution  at  precisely  critical  frequencies. 
Svich  a  solution  may  well  differ  from  the  limit  of  the  non-critical  solution 
which,  it  seems  to  us,  is  the  observable  quantity  here .J 

Jirst  we  shall  discviss  the  limiting  value  of  the  solution  given  by 
(5)  and  (7). 

To  do  so  we  note  that  if  v  is  such  that  sin  0  tends  to  zero,  then  the 
denominator  of  the  v-th  term  of  (5)  goes  to  zero.  The  corresponding  exponential 
factor  in  this  terra  becomes  exp(+  ikx) .  This  analytical  phenomenon  is  well 
knoim  since  the  time  of  Rayleigh,  who  pointed  out  the  connection  between  grating 
resonances  and  the  case  when  a  spectral  order  passes  off  tangentially.  Of  course. 
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the  approximation  f(G)  ^  f^C©)  breaks  down  completely  in  this  care.  However, 
we  shall  show  that  the  use  of  (?)  leads  to  an  intelligible  result.  Our  task 
will  be  to  compute  the  limit  of  f(e  )/sin  ©  as  ©  -^  0,  n.  In  addition  it 
will  be  noted  that  all  other  terms  of  the  series  (5)  must  be  re-examined  in 
this  case,  since  one  or  both  of  the  quantities  A,  and  A_  diverge,  and  the 
quantities  A-fCo)  and  A_f(n)  become  indeterminate. 

In  order  to  proceed  it  is  necessary  to  scrutinize  the  manner  of  di- 
vergence of  A-  and  Ap»  This  divergence  can  occur  in  several  ways  (cf.  (8)). 
(a) :  kd(l-cos  ^  may  approach  an  integral  multiple  of  2n,  or  (b) :  kd(l+cos  JiO 
may  approach  such  a  value,  or  (c):  both  case  (a)  and  case  (b)  may  occiir  siraul- 
taneously.  (For  the  last  mentioned  case  to  occur  tan  (0/2)  must  necessarily 
be  rational.)  Suppose  therefore  that  there  exist  positive  integers  p  ,  p. 
such  that  either 

(11a)         kd(l-cos  0)  =  2p  n  +  5  6  -»  0  , 

or 

(lib)        kd(l  +  cos  0)  -  2p^n  +  6^  ,  6^^-^  0 

or  both.  In  case  (a)  the  denominator  of  the  term  v  »  p  vanishes  in  iS)   and 

cL 

also  A^  diverges.  In  case  (b)  it  is  the  term  v  »  -  p.  which  is  critical  while 
A_  diverges »  It  can  be  shown  by  comparison  with  the  binomial  expansion  of 
(l-z)"  '  that  for  small  6^  .  we  have 


(12)  4         -^      oMovAF.       y^/^ 


a,b 

V-^      2/(2kd6^^^)-/^^ 

where  the  root  is  positive  or  positive  imaginary.     In  case  (c)   (ll)   and  (12) 

show  that 

(13)  liJn  Aj^/Ajj  =     lim    /S^ 


a      • 


■it 

See  Appendix  I-A, 
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If  we  keep  ^  fixed  and  let  k  approach  a  critical  value,  then  the  limit  is  cot  ^/2. 
If  k  is  a  suitably  fixed  value  and  ^  approaches  a  correspondingly  critical  Vcilue, 
the  liniit  is  +  i.  On  the  other  hand,  (6)  allows  us  to  express  k  sin  ©  approxima- 
tely in  terms  of  5  and  6.  ,  when  v  ■  p  and  v  =  -  p,  respectively.  Introducing 
k^=  k  sin  ^  .   we  find 

(II4)     k^  -»  (2kd6^)^^  /  dj      k  ^  ->  (2kd6.  )^/2  /  d 


Pa        ^    '   '       -Pb 


so  that 


(15) 


2 


lim  k   At  »  lim  k    A_  -  t-  , 
Po  ^        -Pv  ^    ^ 


Since  A-  alone  diverges  in  case  (a),  we  can  calculate  the  limit  A^fCO) 

directly  from  (9),  and  then  (7)  becomes 

(16)  f(e)  -  fj(e)  *  J^H^n^  fi(9'  -  A^a3-\(o)  %(«)  • 

We  noticed  earlier  that  for  large  and  non~critical  kd  the  single  scattering  hy- 
pothesis f(e)  "  fjCO)  supplies  a  first  approximation  to  (?)•  This  is  no  longer 
the  case  as  a  critical  value  is  approached.  Indeed,  neglecting  the  correction 
supplied  by  terms  of  order  A_=  (1/  )^  )  in  (16)  we  find 

(17)  f(e)  ->  fj(e)  -  fj^   fj^(e)  , 

which  differs  from  f_(©)  by  a  finite  amount.  In  case  (b)  a  similar  result  is  ob- 
tained. In  (17)  it  should  be  noted  that  the  denominator  is  the  forward  scattering 


See  Appendix  I-B  for  the  evaluation  of  lim(A^/A_)  in  the  various  cases;  also 
for  the  evaluation  of  the  limit",  in  equations  (lU)  and  (15) . 


-  9  - 


ainplitxide  of  an  isolated  cylinder.  Its  negative  real  part  gives  the  scattering 
cross  section,  by  a  well-known  theorem.  The  vanishing  of  the  cross  section  im- 
plies that  2n 


If^C©)!^  de   2   0  . 


This  in  turn  implies  frCO)  ^  0  for  all  6,  whence  also  the  vanishing  (everywhere) 
of -the  field  scattered  by  the  cylinder  when  excited  by  the  incident  field  exp(ikx). 
The  incident  plus  scattered  field  could  not  then  combine  so  as  to  fulfill  homo- 
geneous boundary  conditions.  Consequently  frCo)  cannot  vanish. 

We  now  consider  case  (c).  In  this  case  we  have  to  consider  two  possibili- 
ties: either  a_  vanishes,  or  it  does  not.  Both  possibilities  could  actually  occ\xr. 
It  is  easily  verified,  by  calculating  the  cross  section  explicitly,  that  for  the 
case  of  a  circular  cylinder  a.  cannot  vanish.  Cn  the  other  hand,  for  a  strip  of 
zero  thickness  along  the  y-axis  the  symmetry  of  the  scattered  field  shows  that 
a»  does  vanish. 

It  is  of  interest  to  discuss  the  determinants  cl.  ,  a^  and  a..  These  are 
properties  of  a  single  cylinder.  It  can  be  shown  that  when  a-  vanishes  so  do  ou, 
and  Og.   To  see  this  define  v  =  frCO)  F_(x,y)  -  fjj(O)  ¥Ax,y),   which  is  an  out- 
going wave  function.  Then  form  the  quantity  E(v)  5  Im  d  v(3v/dn)ds  over  the 
body*  Then,  using  the  boundary  condition   for  F,  and  F^,  and  the  well-known 
expression  for  the  far  field  amplitude  as  a  surface  integral,  we  find 
E[v3  =  ijRe  ) 7,(0) a- 1  ,  Consequently  the  vanishing  of  a-  implies  E  [v]  5  o.  Since 
V  is  an  outgoing  wave  function,  it  now  follows  that  v  is  identically  zero.  (The 
argument  is  that  used  in  the  proof  of  the  uniqueness  theorem  for  diffracted  fields, 
which  actually  reqtiires  only  the  vanishing  of  the  energy  flux  E(v)  over  the  surface.) 

See  Appendix  H. 
■^  Thus,  for  E 1 1,  Fj^-  -  e"^^  on  the  body. 
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Having  now  shown  that  a-  s  0  implies  v  2  o,  consideration  of  the  far  field  form 
of  V  shows  that 

fj^(0)fj^(e)  -  fj^(0)fj^(9)  =  0 

for  all  e»  Now  by  the  reciprocity  theorem  for  plane  waves,  we  easily  see  that 
fj(0)  =  f^i^  +  n),   fj(n)  =  fj{^  +  n),   and  fj^(O)  -  fj^(n),  (Recall  for  example  that 
f_(0)  is  the  far  field  amplitvide  in  the  direction  0  produced  by  a  plane  wave  in- 
cident in  direction  ^  while  f_(^  +  n)  is  the  far  field  amplitude  resulting  in  ths 
direction  0+",  when  a  plane  wave  is  incident  from  the  direction  6  =  0.)  After 
these  substitutions  it  follows  directly  that  a. «  0.  The  vanishing  of  ol.  implies 

(by  direct  substitution  of  iL(n)f_(0)/fjj(0)  for  fj(n))  that 

fj(0)a^ 

Consequently  a.  s  0  also  implies  0-=  0,  as  was  to  be  shown,   [it  seems  worthwhile, 
in  passing,  to  inquire  as  to  when  a.  can  vanish.  The  necessary  and  sufficient  con- 
dition (for  perfectly  conducting  bodies)  in  the  case  of  E||  is  that  the  boundary 
is  composed  of  one  or  more  vertical  zero -thickness  strips.  In  the  H||  polariza- 
tion any  number  of  horizontal  strips  may  be  added.  The  spacing  of  the  verticeil 
strips  must  be  according  to  half-integral  multiples  of  the  wavelength.  Thus  in 
the  H 1 1  case  we  can  have  a  rectangle  of  suitable  width.  The  back  and  forward  scat- 
tered amplitudes  agree  for  E] |,  and  differ  only  in  sign  for  H| |.  The  result 
generalizes  well-known  syiranetry  properties  of  a  screen  of  zero  thickness.  If  a 
general  homogeneous  boundary  condition  is  to  be  satisfied,  one  can  derive  a  dif- 
ferential equation  relating  the  eqviation  of  the  curve  and  the  coefficients  in 
the  bo\mdary  condition*  These  results  all  follow  easily  by  applying  the  appropriate 
boundary  condition.  For  the  E||  case,  for  example,  we  have 
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on  the  boundary.  If  a_=  0,  then  v  vanishes  identically  and  hence  x  must  be  con- 
stant on  the  boundary.  Taking  the  origin  of  x  at  a  boundary  point,  we  see 
f^(0)  =  1L(0).  Hence  at  other  boundary  points  sin  kx  =  0,  whence  the  spacing 
can  be  deduced.  The  sufficiency  is  shown  by  defining 

^  =  ^L-  ^R  * 
say,  in  the  case  E| |«  Then  W  will  be  seen  to  vanish  on  the  boundary  and  hence 
W  £  0,  so  that  F.=  Fr,.  Consequently  frC©)  «  %(®^  ^°^  ^^  ®»  whence  it  is  easily 
shown  that  ol.  ,  a^  and  a.  vanish.J 

Returning  now  to  case  (c)  we  first  note  the  result  when  a.  /  0.  Taking 
the  limit  in  (7),  with  f(0)  and  f(n)  given  by  (9),  as  both  Aj^  and  A^  approach  oo, 
we  obtain 

(18)    f (0)  =  fj(©)  -  f^f^Ce)  +  ^2^^^^  /  "3  • 

On  the  other  hand,  if  a^   vanishes,  then  so  do  a-,   and  a„j  and,  using 
the  relation  between  fj(©)  and  £□(©)  we  obtain  in  general 


(19)    f(e)  =  fj(e) 


and  therefore  in  the  limit 


-  [V  \]^l(°> 


(20)    f(e)  =  fj(e)  -  fj^   fj^(e)  . 

We  have  now  calculated  the  limit  which  f(©)  approaches  in  cases (a),  (b) 
and  (c).  This  allows  immediate  calculation  of  the  complex  anqjlitudes  of  the  non- 
critical  spectra,  i.e.,  the  terms  in  (5)  for  which  the  denominator  does  not  tend 
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to  zero.  We  still  have  to  examine  tlie  terms  v  ■  p  ,  and/or  v  =  -  p,  ,  where  p 

and  p.  are  given  by  (ll).  s^ 

In  case  (a)  we  find  ,  using  (l5), 

2f(©  )        ,  ^  ^1 

where  -p(0)  =  a  fj(9)  -  o-i^tC©)  -  S^R^®)  (note  that  p(0)  =  O) .  ^ince  kd  is  always 
large,  this  simplifies  to 

2        '^V  ^l(°> 

(22)  liniH    sin  6  -^     "  fTToT       * 

Pa         ^ 

The  result  (21)  furnishes  a  correction  to  this  formula. 

Thus  the  total  coefficient  in  the  last  spectrum  tends  to  a  finite 
quantity  as  this  spectrum  tends  to  the  horizontal. 

In  case  (b)  we  find,  for  the  cidtical  term  v  =  -  p,  ,  the  corresponding 
approximate  value  -  fj(n)/f  (n)  =  -  fj(n)/f-(0). 

Turning  now  to  case  (c),  consider  first  the  general  situation,  i.e., 

a-/  0.  We  find  the  coefficient  of  the  term  V  =  p_  tends  to 
3  a 

2f(©  )  , 

(23)  l^kd  sine      -^-H"Zr-   • 

^a       -^       -^ 

Similarly,  for  v  »  -  p. 

2   ^(V   _    "2   2  p'(n) 


lim 


E3  sin  9..     ~  a,  ~  kd  a^ 


» 


V     "3   '^^  "3 
which  is  again  a  finite  result. 

As  resonance  is  approached,  in  the  case  a-.=  a^=   <l  =  0  we  get  instead. 


^°^  ^  =  Pa                           2      ^(«v) 

fi(o) 

^^  kd    sin  e^ 

fj^(o)  l+cr| 

See  Appendix  III. 
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plTis  a  correction  terra  in  (kd)~  •     Here   O"  is  the  lindt  approached  by  the  ratio 


Ap/ft_     (cf •   (13)   and  the  subsequent  text) . 


* 


U.  Remarks  on  energy  conservation 

For  any  wave  function  that  is  regular  in  the  exterior  of  the  grating 
elements,  the  total  energy  passing  tlxrough  the  boundaries  of  a  typical  cell  will 
be  equal  to  the  energy  passing  through  the  surface  of  the  typical  cylinder  con- 
tained in  that  cell,  as  can  be  seen  from  direct  applications  of  Green's  theorem. 
(F\xrthermore,  due  to  the  periodicity,  the  energy  passing  through  the  vertical 
boundaries  of  the  cell  is  zero.) 

Consequently,  of  course,  energy  is  conserved  in  the  correct  solution 
just  as  well  as  in  the  incident  field.  For  an  incorrect  solution,  however,  the 
degree  of  conservation  of  energy  will  be  measured  by  the  value  of  the  surface 
integral  of  Poynting's  vector,  provided  we  are  dealing  with  a  suitably  periodic 
wave  function,  which  is  regular  in  the  exterior  of  the  cylinders  of  the  grating. 
Now,  with  regard  to  the  single  scattering  solution,  according  to  which  each  cy- 
linder contributes  to  the  field  as  if  in  isolation,  it  is  easy  to  show  that  the 
boundary  values  are  adopted  vdth  an  error  of  order  A..,  say.  Consequently  energy 
is  conserved  to  order  zero,  except  in  the  resonance  case.  In  that  case  A,  di- 
verges, so  that  energy  is  not  conserved  at  all.   [To  see  this,  we  simply  note,  for 
example,  that  a  cylinder  at  a  distance  nd  to  the  right  of  a  given  "  central"  cy- 
linder contributes  a  field  of  the  order  of 

H  (nkd)  f^-Cn) 
o      1 

to  the  boundary  values  on  the  cylinder  selected  as  "central".]  On  the  other 

hand,  in  our  analysis,  the  contribution  of  a  single  cylinder  is  given  by  formula 

(7a) »  It  is  easy  to  show  that  then  the  stirface  energy  integral  vanishes  to  order 

See  Appendix  IV. 
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(kd)  '     and  that  the  error  is  not  only  of  higher  order  in  (kd)~  for  a  non- 
critical  case,  but  it  is  also  uniformly  bounded,  i.e.,  does  not  diverge  as 
resonance  is  approached.  Consequently  in  our  result  energy  is  conserved  to  the 
proper  order,  and  this  independently  of  the  "  resonance". 

It  should  perhaps  be  added,  however,  that  conservation  of  energy  in 
itself  does  not  seem  to  guarantee  that  the  boundary  conditions  are  even  approxi- 
inately  fulfilled.  From  the  vanishing  of  Im    u  3u/dn  ds  over  a  surface,  only 
the  orthogonality  of  the  function  and  its  normal  derivative  follows.  For  exan^jle, 
the  exact  solution  for  a  grating  of  perfect  dielectric  cylinders  is  certainly 
useless  as  a  substitute  for  the  solution  of  the  corresponding  problem  of  conduct- 
ing cylinders.  Tet  it  certainly  fulfills  conservation  of  energy.  Consequently, 
while  the  energy  law  is  useful  as  a  check,  proposed  approximate  solutions  always 
require  additional  support. 

5.  Comparison  with  other  theories 

Since,  as  is  well  known,  the  theory  of  a  grating  of  conducting  cylinders 
implies  a  corresponding  theory  for  the  conducting  reflection  grating,  it  seems 
useful  to  compare  the  results  so  obtainable  from  our  solution  with  results  ob- 
tained by  Artmann.  Artmann's  results  are  obtained  under  the  assumption  of  pro- 
tuberances small  ccsnpared  to  wavelength,  xfithout  restriction  as  to  spacing.  By 
specializing  his  results  to  the  case  of  large  spacing,  and  specializing  o\rr  results 
to  the  case  of  thin  cylinders,  a  valid  comparison  is  therefore  possible.  However, 
we  require  explicit  expressions  for  the  quantities  fj(0),  etc.,  since  Artmann's  re- 
sults are  expressed  in  terras  of  the  geometry  of  the  surface.  Tq  this  end,  we 
have  calcvilated,  by  the  method  of  small  perturbations,  the  pattern  produced  by 
diffraction  by  a  thin  strip  of  thickness  y  ■  26f (x)  where  ke  is  small  • 

See  Appendix  V  for  this  calculation. 
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We  consider  a  symmetrical  strip  y  =  +  ef(x),  for  |x|  <  X  •  The  incident 
magnetic  field  has  a  single  component,  which  is  perpendicular  to  the  xy-plane  and 
is  given  by  F.  =  exp(ikx) .  Then  if  we  write 

H(x,y)  =  exp(ikx)  +  Fj^(x,y)  , 

Fj^(x,y)  -  eF[^^x,y)  +  e^[2^(x,y)  +  ... 


(23) 
we  find 


.(1) 


(2U)    Y)^'U,y)     -  f 


ikx«  t 


o  o 


(1) 


yCT) 


O'^^i^O  > 


-k 


F[2)(x,y)  --I 


-J? 


(25) 


ay' 


,(1) 


.  H-lk/(x.X^)2.y2Uc^   . 


These  resvilts  are  obtained  by  expanding  the  boundary  condition  for  F- 
in  a  power  series  in  s.  It  is  easily  verified  that  the  function 
exp(ikx)  +  eF^  +  c  ^i   fulfills  the  boundary  condition  to  order  e  •  From 
(2I4)  and  (25)  we  find 

A 


(26) 


fi.(e) 


ke 
T 


-i 


t,     ,  ikz„  (1-cos 
f  (x^)  e   ° 


dx  +  0  i 


(ke)' 


The  leading  term  vanishes  for  6=0,  since  f (+  O  ■  0«  This  is  natural  since 
Re  >f, (0)?   is  quadratic  in  e  by  the  forward  amplitude  theorem  .  To  obtain 
f,(0)  we  use  (25)  together  with  (2I4).  After  various  integrations  by  parts  we 


obtain 
(27) 


h^O) 


-k  6 


00 

2  2   ^ 


'00 


|g^(t) 


dt  + 


o[(ke)2]. 


See  discussion  following  equation  (l?) . 
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where  . 

(28)  0(t)  «=     f  (x^)exprikx^(l-t)1  dx^ 

Here  we  have  used  the  expression 

(29)  H^l\|x|)  -  i:  '     ""^ 


dt, 


where 


-00  /b^-  1 

y^tr-^    =  -  i  )4  -  t^'  for  |t|  <  1  . 

It  is  easy  to  verify  that  Re  if^(0)|   is  negative,  while  Im  >f-(0)i  is  positive. 
Let  us  now  briefly  discuss  the  solution  of  the  reflection  problem.  It 
is  well  known  that  to  obtain  this  solution  from  the  cylinder  solution  we  have  to 
consider  a  grating  of  cylinders  symmetric  with  respect  to  the  plane  y  =  0,  and 
construct  the  function  u(x,y5^)  +  u(x,y{-0) ,  where  u(x,yj0)  is  the  solution  (for 
the  grating  of  cylinders)  corresponding  to  the  incident  wave  exp|ik(x  cos  {/  +  y  sin  JJOJ^ 
(The  upper  sign  is  taken  when  we  have  H| |,  while  the  lower  sign  refers  to  the  case 
of  E||.)  In  calculating  by  means  of  our  formulas  (7),  (8),  (9),  (lO)  we  therefore 
have  to  consider  the  subscript  I  to  refer  to  both  the  directions  of  incidence  0 
and  -  ^«  Now  if  ^  is  replaced  by  -  J^  then  A,  and  A_  are  unaffected,  while,  be- 
cause of  the   symmetiy,  the  quantity  tSoi)   for  the  direction  of  incidence  -  j2f  is 
equal  to  i'-r(-  6)  as  calculated  for  the  direction  of  incidence  +  ^.  Consequently 
also  fy(0)  and  f-.(Ti)  remain  unaltered  under  this  transformation,  whence  the  same 
is  true  of  ol.  and  a^  of  equation  (10).  Thus  in  calculating  with  equation  (7) 
for -the  direction  of  incidence  -  0,  the  first  term,  f-p(6),  is  replaced  by  f-j.(-6) 
and  the  remaining  terms  are  unaltered.  Consequently,  for  E | | ,  the  necessary  sub- 
traction leads  to  a  simple  value  for  the  total  f(e)  in  the  reflection  problem,  i.e», 

^total^^^   =  ^l(«^  -  ^l(-  «^   • 
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Since  fy(e)  is  an  analytic  function  of  cos  6,  it  is  easily  verified  that  there 
is  no  divergence  of  the  terms  f(e  )/sin  Q    as  6  — ^  0,  in  (5),  so  that  there 
is  no  resonance  phenomenon. 

On  the  other  hand,  in  the  case  of  H] |,  we  have  to  add  the  terms  f(6) 
and  f(-  6).   [it  is  interesting  to  note  in  this  connection  that,  if  the  cylinders 
are  taken  to  be  of  rectangular  cross  section  and  of  unit  or  integral  wavelength 
thickness,  then  it  follows  from  our  discussion  of  the  determinants  a  in  a  pre- 
vious section  that  fj{Q)   +   f-r(-  ©)  «  0  if  we  choose  ^  =  n/2.  In  this  case, 
therefore,  our  result  for  the  reflection  grating  is  the  trivial  result 
exp(iky)  +  exp(-iky),  regardless  of  whether  the  spacing  is  critical.  However 
this  result  is  easily  verified  for  any  collection  of  rectangular  protuberances, 
for  any  spacing  whatever,  if  0  =  n/2  and  the  thickness  is  half  a  wavelength.  A 
protuberance  of  half  a  wavelength  corresponds  to  a  cylinder  of  thickness  equal 
to  a  wavelength.J 

We  continue  our  discussion  of  the  H| |  case  by  specializing  to  the  case 
of  thin  cylinders.  For  strips  of  zero  thickness  the  far  field  amplitude  is  well 
knoim  to  be  odd  in  6.  Consequently  the  term  f  j(©)  +  f-.{-   ©)  will  be  of  order 
ke  where  e  is  the  thickness.  If  we  turn  to  the  next  term,  however,  in  the  case 
of  resonance  [cf.  (17)3>  we  find  a  different  order  of  magnitude.  For  the  re- 
flection grating  this  term  becomes  in  fact 

fjCo) 


-  TTToT  fi^")  '  'l^'  '^] 


Kow  f^(0)  =  ^d(0  ■•■  ")  ^7  reciprocity,  -^br  zero  thickness  this  quantity  is  identi- 
cally zero  and  for  a  thin  strip  it  is  of  order  ke.  The  quantity  f-j.(e)  +  fjC-  ©) 
is  also  of  order  ke.  But  the  denominator  frCO)  is  of  order  (ke)   [cf.  (27)3 • 
Therefore  the  total  amplitude  is  a  finite  nvimber,  no  matter  now  small  ks  is. 
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provided  vie  havB  a  critical  combination  of  incidence  angle  and  wavelength.  This 
is  in  contrast  to  the  situation  for  non-critical  wavelengths  •  Artmann  has  shown 
also  that  the  maximum  amplitude  occurs  at  slightly  larger  wavelengths  than 
critical,  and  that  the  shift  is  independent  of  the  order  of  the  spectrum.  To  see 
this  from  our  result  we  return  to  the  formula  for  f(S  )  at  a  non-critical  wave- 
length. Consulting  (7)  and  (9),  and  neglecting  terms  of  order  A_,  we  find,  for 
a  case  when  A^  is  large  but  not  infinite, 

f  (o)f  (e  ) 

(30)       f(9^)  .  f(-  e^)   -  fj(©^)  .  fj(-  9^)  *     (i)a^)  t   f^(o)   • 


From  (12),  A^ '^  2(kd  5  )"  ^  ,  where,  if  5  is  negative,  the  positive  imaginar7 
root  is  taken.  Since  -Re  ^f^(0)t   is  positive  by  the  forward  amplitude  theorem, 
there  can  be  no  cancellation  between  l/A^  and  fT(0)  if  5  is  positive.  On  the 
other  hand,  if  6  is  negative,  then  l/Ay  is  positive  imaginary.  In  (2?)  we  see 
that  Im  \f^(0)i  is  also  positive.  Thus  the  imaginary  part  of  the  denominator 
will  vanish  at  a  certain  negative  value  of  5  ,  given  by  the  formula 


2  2    °° 


(31) 


y^;;^  ._  kV    r    i^(t)rH^(-t)r  at  . 

'    'a'     2n    I       /—K — I 

A      /t-1 


Consequently  [cf.  (lla)J  a  maxiraim  occxirs  for  kd  slightly  less  than  critical,  in 
general,  since  the  numerator  will  only  vary  slightly  with  such  a  small  variation 
in  k.  Equation  (31)  shows  that  5  is  independent  of  the  order  v,  in  which  the 
maximum  is  being  observed,  and  is  also  independent  of  the  index  p  of  the  order 
which  tends  toward  grazing  .  Furthermore,  it  is  easily  seen  that  if  f  is  the 


Then  the  leading  term  is  f-.(e)  +  f^(-e),  which  tends  to  zero  vri.th  the  thickness, 
vJe  keep  kd  fixed,  and  vary  the  angle  of  incidence  in  order  tc  vary  p  . 
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an5)litude  at  the  critical  value,  while  f  Is  the  value  at  the  maxirauin,  then 
(f  -  f  )/5  is  of  the  order  of  f  /6,  i.e.,  inversely  proportional  to  6,  Conse- 
quently the  maximum  is  very  steep. 

It  is  possible  to  show  that  all  these  results  are  in  agreement  xd.th 
the  restilts  of  Artmann  provided  we  replace  Artmann's  sums  by  integrals,  in  ac- 
cordance with  the  fact  that  kd  is  to  be  assumed  large.  Using  our  notation,  we 
can  write  Artmann's  result  (cf»  his  equation  (3.32))  as 


6 


r  ^a  1 

f(x)exp  ix  -T-  +  k  cos  ©  -k  dx 


a 

Using  the  fact  that  d  is  large  we  can  replace  the  above  sum  by  an  integral,  via  the 

substitution  k  cos  G  =  kt,  2n/kd  =  dt,  -co  <  t  <  oo.  3ince  5  is  small  we  ha^/e 

2n(v-p  )/d  «»  kt  -  k,  while  the  term  5  /d  in  the  exponential  can  be  negelected. 
s  a 

This  results  precisely  in  our  formula  (31).  On  the  other  hand  Artmann's  theory  is 
also  valid  for  small  thickness  of  protuberances,  even  if  the  spacing  is  not  large, 
while  our  results,  although  restricted  to  large  spacing  (kd  »  2n),  do  not  require 
that  the  thickness  is  small  compared  to  wavelength. 

ri2i 

We  now  tiorn  to  a  brief  comparison  with  new  results  of  Twersky*-  -■  which 
appear  to  be  in  disagreement  with  oxirs.  To  summarize  this  disagreement,  we  note 
that  in  [l2]  a  functional  equation  is  solved  for  the  amplitude  in  the  reflection 
grating.  This  functional  equation  is  solved  approximately  by  iteration,  the  justi- 
fication for  retaining  only  the  leading  term  being  the  assiimption  of  large 
spacing  (cf.  [l2],  equations (70)  and  also  page  11).  Then  the  result  is  a  quantity 

?  -    7(9  -n-9  )  which  is  in  t\im  defined  as  the  solution  of  a  new  integral  equa- 
vo      v'    o  t.     -I 

tion  (cf.  [12],  equation  (60),  (56))  ^,^ 


1      -     f(e  ,n-e  )  -  f(9,n-©)-i. 
vo      V    o        V*    0    2n 


r 


f(9,n-J2O7(0,n-J2f^)d(2f  , 


; 


-V2 


-  20  - 


in  which  the  quanity  f  is  the  restilt  based  on  single  scattering  theory,  and  f  is 
to  be  solved  for,  with  f  as  given  data  .  It  is  shown  that  f  is,  \inder  certain 
circumstances,  purely  imaginary,  and  that  consequently  if  one  retains  terms  of 
higher  order  of  magnitude  in  the  specular  amplitudes,  while  neglecting  them  in 
the  other  spectra,  then  energy  is  conserved  exactly  (  [l2] ,  pp»  29-30).  It  is 
then  en^Dhasized  in  ([l2]],  p.  30),  that  the  results,  (which  according  to  page  29 
and  eqiiation  (70)  are  obtained  for  large  kd)  "  differ  significantly  from  a  single 
scattering  approximation  in  that  they  contain  f  instead  of  f ".  Arguments  are 
given  to  show  that  the  presence  of  f  rather  than  f  would  be  incorrect  physically. 
The  results  insofar  as  they  are  explicit  thus  seem  to  be  in  disagreement  with  oia* 
own  even  before  the  passage  to  the  limit  of  resonance,  since  o\rr  naive  leading 
term  is  the  single  scattering  term,  in  the  non-critical  case,  and  our  principal 
resiilt  is  the  correction  to  this  term,  so  that  in  the  limit  of  infinite  spacing 
our  results  reduce  to  single  scattering.  We  reserve  further  comment  ;intil 
Twersky's  results  have  appeared  in  final  formj  but  it  may  be  noted  that  the  dis- 
crepancy would  be  removed  if  the  validity  of  Twersky's  result  were  additionally 
restricted  to  small  scatterers,  or  at  any  rate  to  scatterers  whose  scattered 
airQslitude  is  small,  for  then  his  integral  equation  shows  that  f  =  f .  It  there- 
fore seems  possible  to  us  that  the  condition  for  the  usefulness  of  Twersky's 
iterative  solution  of  his  functional  equation  (70)  actually  dpends  on  the  smallness  of 
the  single  scattered  amplitude,  rather  than  the  large  kd  he  actually  postulates. 
Then  there  would  be  no  discrepancy  in  the  limit  of  large  spacing.  The  inherently 
plausible  results  proposed  in  {l2],   on  a  "  heuristic"  basis,  for  the  case  of  one 
mode  of  propagation  would  also  be  validated  in  a  logical  sense.  The  reason  is 

In  tMs  formula  the  superscribed  bar  does  not  mean  complex  conjugate. 
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that  for  the  reflection  grating  our  suggestion  would  imply  the  assumption  of 
the  smallness  of  the  single  scattered  amplitudes  of  a  typical  isolated  pro- 
tuberance. Since  this  is  the  case  for  protuberances  of  thickness  small  com- 
pared to  wavelength  the  results  should  then  agree  vdth  those  of  Artmann 
regardless  of  spacing* 
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Appendix  I-A 

Eqiiation  (12)  is  proved  as  follows.    We  have  from  (8),  (11), 


But  the  binomial  forravila     gives 


00     (^)''(-i): 


(i.z)-^/2  -  21 

^0       (-|-r)lr! 


T 


Using  the  identities   ["""(z+l)  «=  z   {~Xz) ,    (~\*  +  z)  [~tx  -  z)  ■  n/cos  n  z,  together 
with  Stirling's  formula 


we  find 


or 


U..,-iA^   i   f    ^ 


.'-['vr-  p^  ?  =1^ 


and 


For  sTifficiently  small  5^  ^: 
and 


\,R  ~    /^ 
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Appendix  I-B 

We  evaluate  lira  A./ZL,  by  vising  the  expressions  (11a),  (lib)  and  (13) • 

Suppose  first  that  0  is  fixed  and  that  k  — >  k  ,  a  given  critical  value.  But 

ftora  (11a),  (lib): 

5^  kd(l  +  cos  0)  -  2pj^n 

37     "     kdtl-cos  ^)  -  2p^n 

As  k  -»  k  ,  the  numerator  and  denominator  of  the  expression  on  the  right  approach 
c 

zero  (this  is  the  definition  of  a  critical  frequency  in  case  (c)).  Using  L'Hopital's 

rule,  we  find 

5     1  +  cos  0 


k  — >  k   a 

0 


cos 


lira   T-  ■  ^ Ai    "  cot 


and 

^  o  li^    yV^  -  cot  I  . 


lim 
k  ->  k 


In  the  case  where  k  is  fixed  (i.e.,  at  some  value  which  is  critical  for 

a  certain  value  of  0)  while  jZf  approaches  the  critical  angle  0  ,  we  again  use 

L»H8pital's  rule  to  obtain 

&.      kd  sin  0 

r^r^    K    '     '   kd  sin/  '  "  ^ 
0  — >  Of   a  '^c 

c 

and 

A 


L 
c  ^    —  ^c 


0  ->  0^  ''r    i^  -»  0« 


Equation  (lit)  is  derived  as  follows.  Substituting  v  «  p  in  (6),  and 


a 


using  (11a),  we  have 


JETB^ 


The  case  v  =  -  p,  is  identical.  Equation  (l5)  follows  from  (12)  and  (lii). 
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Appendix  II 

It  is  required  to  show  that 


Im  d)   V  ^  ds  »  U  Re  ^J^lCo)  a.' 


{■ 


where 

V  =  fj^(O)  Fj^(x,7)  -  i^(0)Fj^(x,y) 

and  B  is  the  boundary  of  a  given  isolated  cylinder. 

We  begin  by  expressing  F(x,y),  F(e)  (see  (3),  (U))  as  surface  integrals. 
The  response,  or  scattered  field,  F(x,y),  is  given  by 


F(x,y)  -  u(x,y)  -  \j^^^U,7) 


'B 

where  r  p  is  the  distance  from  an  arbitrary  point  on  B  to  the  field  point  P(x,y)j 
the  integral  follows  from  applying  Green's  theorem  to  two  reduced  wave  equations 
(see  [2]  for  details).  Then 

Fjj(x,y)  =  \i^,y)     -     e"^^ 

i  rf   ^^^^^  „(1)  rw  ^H 


F 


ikx 


j^(x,y)  =  VL^(x,y)  -  e 

"•iloc   iloc 
where  u^,  u,,  are  the  total  fields  resulting  frou  the  incident  fields  e    ,  e 

respectively.  The  far  field  form  of  F(x,y)  is  then  obtained  by  substituting  for 
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o 

tegral.  Applying  the  definition  (Li)  we  obtain 


H^  (kr  _,)  its  asymptotic  form,  and  suppressing  higher-order  terms  in  the  in- 
o    sP 


so  that 


; 


r      ^.R^^^       -ikp(s)   cos    [©-©(s)] 
D      ^  e 

B 


fr  0(0)     =    r    "^      fe e     "^^   '  L        ^  '-J  ds 


and 


^l,r(°>  '  n  t  -is—  ^      ^^  ' 

.       r  V  i     I  ^L.R^^^       ikx     , 

^l,r(°^  =  IT  j  — ^T-    «         ^^  ' 

\r^"^  =  IT  J  ~^^  • 

The  proof  of  the  required  identity  is  now  a  direct  computation.  B7 

definition 

V  =  fj^(O)  Fj^(x,y)  -  fj^(O)  Fj^(x,y) 


^l(°)  (y  ^"'^)  -  %(°)  (v  -'"")  • 


Since  Upi  'u-r  vanish  on  B,  we  have 


v(s)     =  -fj^(O)   e-^^  +     i^(0)   e^^ 


and 


^  -  V°)  ^  (V  »^'j  -  ?r(°>  It  (v  ^■''°)    . 


so  that 
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/   X  av(s)  |„  /^m2     -ikx     3     /■-        ikx> 


"sr 


SE 


(V 


I-  /„\|2     ikx       8     /—        -ikx\ 


.    f,(0)  1^(0)  e-il-    I    (V  »■'■"' 


.    f,(0)  ?,(0)  e^^      I,    fv  .^^) 


Notice  that  the  sum  of  that  part  of  the  last  two  terras  on  the  right  which  does 
not  include  diL/dn  or  duL/Bn  has  an  imaginary  part  of  Oj  therefore  these  terms 
may  be  omitted  in  forming 

Im 


/B 


In  the  remaining  terms,  those  which  do  not  include  diL/dn  or  3Uj./3n  vanish  upon 
integration  over  B,  It  follows  that 


i 


v(s) 


3v(8) 
dn 


ds  =  -  |fj^(0)r  <P  e 


-ikx 


ikx 


aiu. 


ds 


-  |f^(0)r  t  e^'°'  ^    ds  *  f,(0)1^(0)  t  e-l-^ 


ds 


B 
.     %(0)  7^(0)  I  e^^  ^  ds  • 


'B 


F^ora  the  definitions  given  above,  we  therefore  have 
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^i^«     -    r    j   ^l(0>  V°)   V")  -   \'L^O)fl^in) 


'B 


f,.(o)    J  7,(n)  7„(o)  -  7,(0)  7,( 


r    ^L^^'      S  ^L^"^   ^R 


V^'   "R 


■'} 


h±    f^Co)  Oj 


where  a.  is  given  by  (10)«     The  resvilt  is 


Im^  v^ds  -  URe  J  7j^(0)  a  >   • 


; 


B 


du 


A  sirailar  treatment  can  be  given  for  the  case  of  the  boundary  condition  ^  =  0 
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Appendix  III 

In  case  (a) ,  the  coefficients 

^     f(a  ) 

kd  sin  ©   ^'^^v^ 

may  be  computed  from  (16),  for  v/p»  Ifv=p,  the  denominator 

kd  sin  ©   «  k  d 
Pa    Pa 

approaches  0,  and  the  value  of  the  coefficient  of  index  p  is  apparently  indetermi- 

cl 


nate.  But  from  the  definition  (6)  of  the  spectral  angles  6^^,  we  have 

-1 


f 


©  =  cos 


2«v  _,    ^1 


and  from  the  definition  (lla)  of  p  ,  we  have 

©  =  cos"^  (1)   =  0  . 
Pa 

The  coefficient  of  the  term  of  index  p  is  therefore  given  by 
f(e^  )  f(©^  ) 

Pa         ^      Pa     ^      Pa 
where  the  last  expression  on  the  right  follows  from  (l5)»  Ifeking  use  of  (?)  and 
(9),  and  defining  for  convenience 

D  =  1  -  [Ij^  fj^(O)  .  L^   fj^(n)]  .  A^A^  a^ 

(a.  is  defined  by  (10)),  we  find 


^a 


i  ^1?  h\^  *  4  ¥\>[^i(o)  -  \°i]*\^^(Vt^^^"^"^''2]j 
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Rearranging,  we  have  four  terns  to  consider,  as  6  ->  0: 


,-1 


where 


(i)   D--  A^  fj(e  ) 


(ii)   D-1  a2  L(o)  f^(ep^)  -  f^(0)  fj(Gp  A   , 
(iii)  D-1  A^Aj^  L(n)  f^i^^y-f^W   fjCOp  A   , 


(iv)  -D-^A^Aj^pCe  )  , 


-  P(e^)  =  a3  fj(9^  -  a^  f^(e^)  -  a^   fj^(e^ 


Then  dividing  the  numerators  and  denominators  by  A-  and  letting  ©  — >  0,  we  find 

^a 
for  terms  (i),  (iii),  the  limits: 

(i,a) 


fjCo) 


(iii,a) 


Proceeding  similarly,  we  find  for  (ii)  the  limit 


lim 


fj(o)  fj^(e  )  .f^(o)  fi(e  )K^ 


\  ^3  -  ^l(°> 


which  may  be  rewritten  as  a  product  of  limits : 

1 


j^H^TT^  ,":;„(v^"V«'i 


e .^^^  ii5-e-  • 

Pa  ^*  Pa       ^* 
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But  the  first  lindt  in  the  prodiict  follows  from  (1$),   the  second  follows  from 

L'Hopital's  rxile,  and  the  third  approaches  1.     The  result  is  that   (ii)   has  the 

limit 

2        fj(0)  f^(0)   -  fj^(O)  fj(0) 
^^^'^^  H      k^a^  -     fj^(o) • 

Term  (iv)  is  evaluated  in  the  same  way:  using  p(0)  =0,  the  limit  is 

t 


(0) 

^^^'^^    "k   Vr  ^L^^^    * 

Adding  (i,a)  -  (iv,a),  we  obtain  the  stated  result  (21)  for  the  limiting  value 
of  the  critical  term  of  the  series  (5).  The  corresponding  analysis  in  case  (b) , 
where  A_-^oo  and  the  critical  term  has  index  v  =  -  p,  ,  is  almost  identical;  and 
only  slight  modifications  are  required  to  treat  case  (c),  where  both  A,  and  A^ 
diverge  and  there  are  two  critical  terms. 
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Appendix  IV 

A  -typical  cell  of  the  grating  is,  say,  the  x-interval  (a,a+d)  which 
contains  the  zero-th  element:  the  boundary  B  of  the  element  intersects  the 
x-axLs  at  X  »  a,  X  =  p,  where  a<a<0<p<  a+d.  The  total  energy  passing 
through  the  vertical  boxmdaries  of  the  cell  is  given  by  the  integral  of  the 
Poynting  vector  P,   P  ds,  taken  over  the  lines  x  =  a,  x  =  a+d.  Substituting 
for  P  =  Ex  H  its  time-average  P-  =  ■^^  Re(E-*  H),  and  calculating  H  in  terms  of 
E  fi-om  the  relevant  Maxwell  equation  for  the  time -periodic  field,  we  have  as 
a  measure  of  the  flow  of  energy 

Im  J  u  §  ds, 

where  the  integral  is  taken  over  the  two  lines  x  =  a,  x  =  a+d.  But  from  the 

periodicity  of  the  field  (see  (3)),  we  have 

/  ^ ,  \      ikd  cos  Q^  /   X 
u(x+d,y)   =  e       ^   u(x,y} 

and  taking  normal  derivative  in  the  outward  direction,  we  obtain 

8u(x+djy)   _^  _  -ikd  cos  J2f  du(x,y) 
dx       ~  dx 

It  follows  that  the  integral  over  the  vertical  boundaries  of  a  typical  cell 
vanishes:  /■ 


3u  ,      - 

u  ^  ds  =  0  ♦ 


Im 
J 
But  the  field  exterior  to  the  element  B  is  regular  throughtout  the  cell; 


it  follows  from  Green's  theorem  that 

Ijm    u^ds  =  Im(|)u^ds 

where  C  consists  of  the  vertical  boundaries  x  =  a,  x  =  a+d,  plus  the  horizontal 
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boundary  segments   (at  y  =  +00)   of  the  cell.     Since  the  integrail  over  the  verti- 
cal boundaries  vanishes,  we  have 

3u     ,  T        I  du      , 

u    5JJ     ds     -     Im  u    ^    ds 

where  the  integral  over  H  (the  horizontal  boundary  segments)  measuresthe  total 
energy  flux  through  the  cell.  Now  the  integral  over  C  or,  equivalently,  the 
integral  over  H,  nnist  vanish  for  an  energy-conserving  solution,  and  therefore 
the  integral  over  B  must  vanish  for  am  energy-conserving  solution.  This  is 
trivially  the  case  for  the  correct  solution  (which  vanishes  on  B) ,  but  may  not 
be  for  a  given  approximate  solution.  An  acceptable  approximation  to  u  should 
conserve  energy  to  ttithin  terras  of  an  appropriate  order  (although  this  alone 
does  not  validate  the  solution),  while  a  proposed  approximate  u  such  that 
Ira  d)  u  ^  ds  is,  say,  \mbounded  is  obviously  inacceptable. 

For  the  single -scattering  approximation,  i.e.,  the  assumption  that 

F(x,y)  =  Fj(x,y) 

(cf.,  however,  our  approximation  (7^)),  we  discuss  the  energy  balance  as  follows. 
The  field  on  B  consists  of  the  sum  of  four  fields: 

(i)    the  incident  field  vl-j 

(ii)   the  response  field  F,  (which  cancels  the  boundary  value  of  u^ 

exactly)  j 
00 
(iii)  a  sum  ^^    F_(x+nd,y)  of  response  fields  from  cylinders  to  the 
n=l 
right  of  B| 
00 
(iv)   a  sum  ^    F-j.(x-nd,y)  of  response  fields  from  cylinders  to  the 

left  of  B. 
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The  net  field  on  B  is  therefore  the  sum  of  fields  (iii)  and  (iv),  i.e., 
of  the  approximate  fonn  (see  [2],  4  3,  for  details) 


Forming 


-  A^  fj(n)  e^^  -  Aj^  fj(0)  e"^^ 


Im  (P  w  T--  ds 

'B 


where  w  is  the  field  obtained  by  the  single- scattering  approximation,  we  see 

-1/2 
that  energy  is  conserved  to  order  A^  .  =   0 (kd  '  ) ,  at  non-resonance  fre- 
quencies but  is  not  conserved  at  resonance  frequencies,  since 

w  -  -  Aj^  fj(n)  e^^  -  A^  fj(o)  e"^^ 

on  B  and  either  A^^  or  A.,  or  both  diverge  at  resonance. 

In  the  present  approxim?.tion,  the  field  on  B  is  again  the  sum  of  four 
fields,  as  follows: 

(i)    the  incident  field  u_j 

(ii)   the  response  field  F(x,y)  (cf.  (7a))} 

(iii)  a  sum  J^    H   (kr  )g(r  ,6  ),  i.e.,  of  cylindrical  waves  of 


1 

amplitudes  g(r  ,0)5 
n'  n 


f^iQ  )    f.(e  ) 

^  n     r         2 

^       ^n       - 


Here 


(r  ,©  )  are  coordinates  of  the  field  point,  measured  from  the 
n  n 

origin  (x+nd,G)  of  response  fields  from  cylinders  to  the  right 

of  B  (in  terms  of  OTir  previous  notation  we  are  representing 

r(x-nd,y)  as  an  asymptotic  series  in  decreasing  powers  of  r  ); 

ri  (1) 
(iv)   a  sum  ^/J  H^  '(kr  )g(r  ,©  )  of  response  fields  frcan  cylinders 

-00 
to  the  left  of  B, 
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•t'he  analysis  leading  to  (7a)  involved  replacing  H   (kr  )  by  the  first 
term  in  its  asymptotic  expansion,  and  replacing  g(r  ,©  )  by  an  amplitude  f (G  ) 
independent  of  r  : 

(A.l)        H^l>(la.^)  ^     /it:  exp[i(kr^-  nA)]  , 

(A.2)         g(r  yQ  )  ^     f(e  )  . 
&\  n»  n      ^  n 

This  amounts  to  replacing  the  sum  of  fields  (iii),  (iv)  by 

k^   f(n)  e^^  +  k^  f(0)  e'^^  . 

Adding  (i) ,  we  see  that  the  sum  of  fields  (i),  (iii),  (iv),  is  exactly  cancelled 
by  (ii)}  i«e.,  the  net  field  on  B  vanishes,  to  id. thin  the  degree  of  approximation 
represented  by  (A.l)  and  (A.2), 

Introducing  higher-order  terms  in  the  expansions  of  fields  (iii),  (iv) 
will  of  course  lead  to  a  non-vanishing  field  on  Bj  it  remains  to  be  shown  that 
energy  fas  measured  by  Im  ())  u  (3u/dn  )dsj  is  conserved  to  within  o  (  (kd)""'  j. 
We  accordingly  introduce  the  full  asymptotic  expansion  of  H   (kr  )  in  place  of 
(A.l),  and  an  expansion  in  powers  of  (r  )   in  place  of  (A.2),  viz: 


(B.1)         H^^^kr^)     -     y^exp[i(kr^.  nA] 


where 


00       (-1)\ 

1  +  21   1 

k=L     (2ir  r 


» 


f  (e)        f  (©) 
(B.2)  s(r  ,e  )  -- f(e)  +4 —   ■*■    "^—5   *    •••    • 

(r^) 

The  fvinctions  f-,(e),  f2(©),  •••  can  be  deduced  from  f(e)j  and  if  we  define 
functions  f^,{9),   f,  ,(©),  fL  ,(©),  ...  it  follows  that  each  of  these  has  the 

i,i.       LfX  K,± 
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same  form  (cf.  (7))  as  f(e).  Using  (B.2)  instead  of  (A.2),  attending  for  the 
moment  to  the  second  term  of  (B,2)  and  approximating  r  by  nkd,  we  obtain  for 
(iv)  an  expansion  whose  leading  term  is  a  sum  of  terms  of  type  (nkd)  '  ,  e.g. 

00    inkd(l-  cos  ^ 


fi(o) 


?     (nkd)^/'^ 


This  seizes  converges  absolutely,  even  at  critical  frequencies,  and  an  argument 
of  the  type  given  in  Appendix  III  shows  that  f, (O)  is  finite  at  critical  fre- 
quencies. Therefore  the  use  of  (B,2)  in  place  of  (A.2)  yields  conservation  of 
energy  to  within  o[(kd)"-^/^J  . 

It  follows  that  W3  need  only  consider  further  terms  of  2~!i  ^   (kr  )f(e  )• 

Applying  (B.l),  and  using 

2 
^n  ^  ^^  -^  nd)  *    ^ 

^'n  -~   nd 


f(e^)  -   f(0)  +  ^  f '(0) 


-3/2 
we  again  obtain  a  series  which  is  dominated  by  ^^^  (nkd)  '  ,  multiplied  by  terras 

which  remain  finite  at  critical  frequencies. 


-  36  - 


Appendix  V 

The  pertiirbation  calciilation  leading  to  (2li),  (2^,  is  as  follows. 
According  to  (23),  the  boundary  condition  3H/dn  =  0  on  the  magnetic  field  has 
the  form 

where  the  normal  derivatives  are  evaluated  on  y  =  +  6  f(x)j  or  since  the  symmetry 
of  the  problem  permits  vcs   to  work  in  the  half -space  y  >  0,  we  ;vill  take  the  normal 
derivatives  on  y  =  +  e  f(x).  Then 

on  |dn|  y  ds  x  ds 


■f77\7IaY      )4  *  e^[f'{x)]^ 


The  boundaiy  curve  y  =  e  f (x)  is  assumed  to  satisfy  two  conditions: 

(a)  f  (O)  =0  at  the  leading  and  trailing  edges,  i.e.,  the  curve  fits 
smoothly  into  the  x-axis  (under  this  condition,  there  are  no  singxilarities  in 
the  analysis )j 

(b)  the  expansion 

is  valid,  i.e»,  f  (x)  is  finite. 


3 

order  6  ) 


Next,  evaluating  H  ,  H  ,  on  the  boundary,  we  find  (to  within  terms  of 
X   y 
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an 

- 

dH(x,e  f(x)) 
ay 

= 

aF^^^(x,o) 

'         3y 

3H(x,ef(x)) 
3x 

= 

ike^^.     e 

+   6 


2  ,  aF[2)(x,0) 


ay 


+  f(x) 


a2Fj^^x,o) 


ay 


dx  dx. 


Using  condition  (b)  above,  we  find 


H 


^ 


T7\PZo7 


^ d^ *     ^    ^ 


aF[2)(x,o) 


+     f(x) 


3^j^(x,0) 


Hjj.6  f   (x) 


)f 


TTFS^ 


t                 ikx           2     t         3F^^^(x*0) 
e  f  (x)  ik  e^^  +     e^  f  (x)  — i-^ 


The  boTondary  condition  8H/8n  ■  0  therefore  yields  the  two  conditions  (obtained 
by  setting  the  coefficients  of  e,  e  ,  respectively,  equal  to  O) : 

,(1), 


-^ f  (x)  (ik  e^^J 

aS'[^^(x,o) 


aF^^^x.o) 
ay 


+  f(x) 


ay'' 


.  f  (X)  -^  -  0  . 


Since  the  quantitites  F^  (x,y),  Fi  (x,y)  must  satisfy  the  reduced  wave  equation, 
they  are  determined  as  the  solutions  of  boundary- value  problems  by  these  two  con- 
ditions* 

The  boundary- vaxue  problem  for  Fi  (x,y)  is: 
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ax 


Qy 


a  F 


(1) 


^  -  ikei'^f'(x),    |xl  -i^ 


0  , 


>  i   . 


Applying  Green's  theorem  over  7  "  s  f (x) ,  we  have 


r 


L  an 


jr^^i^ZiG-GAF^iHds  -  (t  M'^g-G^^ 


aF 
"ST 


(1) 


ds 


(where  G  is  the  appropriate  Green's  function),   or 


F|^^\x,y)      "     -  ik 


A 


ikx 


f  (x  )  G  ds     • 


; 


'4 


The  appropriate  Green's  function  (for  y  ^  0),  whose  normal  derivative  vanishes  on 
the  strip,  is 

G 


.  I  h("    (.  /(x-v^^   iyyf)    . 


Substituting  the  value  of  G  at  y  =  0,  we  verify  (2U) ; 


F^^>(x,y)  .  I 


n 


'i 


e   Of  x^)  H 
o   o 


^^^  ^^  /(x-x^)^+  y2'  I  dx^  . 


Solving  the  boundary- value  problem 


ax 


ay 


3>P^(x,o)  ,       aF[i^ 

L  _     ,     f'(x)   -^i f(x)   li 


a2Fj^\x,o) 


ay 


ax 


sy^ 


0    , 


|x|>i? 
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in  the  same  way,  we  obtain  (2$): 

L 


F[^^(x,y)  -  -  2 


^^   3F,^2^(x,0) 


1 


4 1 


3y       o 


aF[i)(x,o)      ^^i^^(^o^ 


dy^ 


d„a)(.yc;T?7)-„. 


To  obtain  (26),  we  recall  (cf.  (U))  that  fj^C©)  is  the  coefficient  of 

yr"^  i(kr-iiA) 

n 
in  the  far  field  form  of  F^  (x,y)  »  6  F£''-^(x,y)  +  eT£  (xjy).  Using 


where 


^  2^ 

XX        X 

r     2r 


(r  cos  e)x 


—  r  - 


-  r  -  X  cos  ©  , 
o       ' 

—1  /2     —1  /2 
and  approximating  p  '  by  r  ■  ,  we  have  from  (2U),  (25): 


f^(e)  - 


k6 


expjikx  (1  -  cos  6)^  f  (x  )dx  +   o((ks)  )  • 
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